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We develop a methodology that enables us to extract accurately the electromagnetic A form factors 
and their momentum dependence. We test our approach in the quenched approximation as a 
preparation for a study using dynamical fermions. Our calculation of the four form factors covers 
pion masses between about 410 MeV and 560 MeV on lattices with a size of 2.9 fm and a lattice 
spacing ^ 0.09 fm. From the form factors we are able to obtain estimates of the magnetic moment 
and the charge radius of the A, which we compare to existing experimental and theoretical results. 
Our non-zero result for the electric quadrupole form factor signals a deformation of the A, pointing 
to an oblate charge distribution. 
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1. Introduction 

Electromagnetic form factors probe the structure of hadrons, yielding information on their 
size, shape and magnetization. While the nucleon form factors and the N ^ ^ transition form 
factors have been studied quite thoroughly both experimentally and on the lattice [[I]], much less 
has been done for the A form factors. Experiments are notoriously difficult due to the short mean 
life time of the A of only about 6 x lO^^^s. Nevertheless the magnetic moments of the A+ and 
A++ 0] have been measured. On the lattice the electromagnetic form factors of the A baryon 
have been evaluated at one value of the momentum transfer in a pioneer quenched study in the early 
nineties [^]. The main advantage of the fixed-sink approach that we adopt here, is that we are able to 
calculate the form factors for all values and directions of the momentum transfer q simultaneously, 
allowing for an increased statistical precision. Furthermore, we construct optimized sources that 
isolate the suppressed form factors enabling us, for the first time, to extract the electric quadrupole 
accurately. Other improvements comprise reduced systematic errors due to finer lattice spacings, 
larger volumes and smaller pion masses. 



2. Lattice techniques 

We are working in Euclidean space-time throughout the whole paper. We use Wilson fermions 
and the standard Wilson plaquette gauge-action. At all the pion masses considered in this work the 
A is a stable particle. 

2.1 Electromagnetic form factors of the A baryon 

The matrix element of the electromagnetic current, , between two A-states can be decom- 
posed in terms of four independent covariant vertex function coefficients, a\{q^), a2{q^), c\{q^) 
and C2{q^), which depend only on the momentum transfer squared q^ = {pf — piY 



{^{Pf,sf)\v^\^{Pi,si)) 



1 Ua(pf,Sf)&^^'u,[pi,Si) 



2mA 



+ 



(2.1) 



and niA denote the energy and the mass of the particle, pj (pf) and Sj (sf) are the initial (final) 
four-momentum and spin-projection, while P = pf + pi. Every vector-component of the Rarita- 
Schwinger spinor Ua{p,s) satisfies the free Dirac equation. Furthermore, the following auxiliary 
conditions are obeyed 

7aM^(/7,5) =0, pau''{p,s)=0. (2.2) 

The vertex function coefficients are linked to the phenomenologically more interesting multipole 
form factors Geo, Ge2, Gmi and Gm3 by a linear relation [|6|]. The dominant form factors are the 
electric charge, Geo, and the magnetic dipole, Gmi, form factors. 
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2.2 Interpolating fields 

We use an interpolating field that lias tlie quantum numbers of the A+ baryon 



(2.3) 



where C is the charge conjugation matrix. To facilitate ground-state dominance we employ a co- 
variant Gaussian smearing [Q] on the quark-fields entering Eq. ( 23 ) 



H{x,y;U) = ^ (t/^(x,0%-A + A,0%+a) 



(2.4) 
(2.5) 



Here q is the local quark field (i.e. either u or d), q is the smeared quark field and is the SU (3)- 
gauge field. For the lattice spacing and pion masses considered in this work, the values a = 4.0 
and n = 50 ensure ground state dominance with the shortest time evolution that could be achieved. 



2.3 Correlation functions 

We specialize to a kinematical setup where the final A-state is at rest {pf = 0) and measure the 
two-point and three-point functions 



Garir\p,tf-ti) = Y^e-'^'fPrl,^{Xoa{tf,Xf)Xra'{U,d)) 

Ga^ir\q,t) = j;^e'''^rl,„{Xaa{tf,Xf)V^'it,x)x,a'{tuO)), 

X,Xf 



(2.6) 
(2.7) 



where is the symmetrized, conserved lattice electromagnetic current. We work with a repre- 
sentation of the Clifford-algebra in which 74 is diagonal. In this representation our choices for the 
F-matrices are 



. 1 / a(*^) \ ^ ^ 1 / 1 \ 

^ =2(^0 oj ^ = 2(^0 oj' 

with k = I,. . . ,3 and a'^^^ being the Pauli matrices. The ratio 

G„\{r,q,t) I Gkk{r\pi,tf-t)Gtk{r\d,t)Gkk{r\d,tf) 



(2.8) 



Gkkir\o,tf) V Gtkir\o,tf-t)Gkkir\pi,t)Gtkir\pi,tf) 



(2.9) 



with implicit summations over the indices k with k = I, . . . ,3, becomes time independent for large 
Euclidean time separations tf — t and f — f,-: 



R^\{r,q,t)^n^,{r,q) 



+ ■ 



niA 



tr 



(2.10) 
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The traces act in spinor-space and the Euclidean Schwinger-Rarita spin sum is given by 

-if + m^ 



2mA 



7ct7t , .PaYz-PxYa 
Oa-c ^ I- 



3m2 



3mA 



(2.11) 



Since we are evaluating the correlator of Eq. ( ^?7| ) using sequential inversions through the 
sink a separate set of inversions is necessary for every choice of vector and Dirac-indices. The 
total of 256 combinations is beyond our computational resources, and hence we concentrate on a 
few carefully chosen combinations given below 



n?® = i:n/,(r4,^) 



(2.12) 
(2.13) 
(2.14) 



From these all the multipole form factors can be optimally extracted. For instance ( |2.12| ) is propor- 
tional to Gm\, while ( |2.14| ) isolates Gei for = 4. 



2.4 Data analysis 



For a given value of the combinations given in Eqs. ( 2.12 ) to ( 2.14 ) are evaluated for all 
different directions of q resulting in the same q^, as well as for all four directions pL of the current. 
This leads to an over-constrained linear system of equations, which is then solved in the least- 
squares sense yielding estimates of Geq, Gei, Gm\ and Gm3- This estimation is embedded into 
a jackknife binning procedure, thus providing statistical errors for the form factors that take all 
correlation and autocorrelation effects into account. 



3. Results 

3.1 Simulation parameters 

The calculation has been performed on a lattice of volume xT = 32^ x 64 in the quenched 
approximation. The lattice spacing has been estimated using the nucleon mass in the chiral limit [|l]] 
to be a = 0.092(3), so the spatial extent of the lattice is ~ 2.9 fm. We work with two degener- 
ate Wilson valence quarks. Isospin symmetry relates results obtained for the A+ to those for the 
A++,a0 and A" since they differ only by a charge-factor. The pion and A masses corresponding 
to the three values of the hopping parameter, k, considered here are summarized in Table [1} We 
use 200 well separated gauge configurations for the calculation. Contributions from disconnected 
diagrams are neglected. The numbers and figures below are obtained from the connected contribu- 
tions to the full electromagnetic current. They also can be interpreted as the iso-vector form factors 
with the iso-vector current = uj^u — dy^d. 
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3.2 Electric charge form factor Geo 

Our results for Geq are displayed in the left panel of Fig. [T[ The momentum dependence of 
our data is described very well by the dipole ansatz 



1 



{l+CEoq^y 

Non-relativisticaUy the slope at = is related to the electric charge radius via 



(3.1) 



(3.2) 



72=0 



for which the results are collected in Table [T[ An extrapolation of the form factor to the physical 
pion mass, linear in m^, is used to obtain the last row of the table. One could consider using 
chiral perturbation theory for the extrapolations This will be done when we have results using 
dynamical fermions. 

After a chiral extrapolation, linear in m^, the authors of ref. obtained 0.63(7) fm, which is 
compatible with our value of 0.691(6) fm. 

3.3 Magnetic dipole form factor Gmi 

Also Gmi {q^) can be fitted to a dipole ansatz 

2n _ aMl(O) 



G 



'-) 



(l+CMi q^Y ■ 



(3.3) 



From the fit parameter qmi the magnetic moment = Gmi (0) 2ml j '•^^ ^ determined. 
It is given in Table |T| Fits and data are shown in Fig. [l[ 
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Figure 1: The two dominant form factors for the three pion masses as well as the result of a linear chiral 
extrapolation to the physical pion mass. The left panel shows the electric charge form factor, Geo, ™d the 
right panel shows the magnetic dipole form factor, Gmi- The lines are fits to the dipole forms Eqs. (3.1) 
and (3.3). The band at = indicates the experimentally measured value as quoted by the particle data 
group. 
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The magnetic moment of the has been measured in several experiments. However not 
all experimental values are in agreement with each other, so the particle data group quotes a rather 
broad band of S.V/Ia? to V.S^Ua? IQ]. The most recent experiment quotes /Xa++ = 6.14(51)/iAr. 
Our estimate of /Ia++ = 2/Ia+ = 4.64(32)/iAr lies in the middle of the error band given by PDG and 
below the result of An experiment that measures the magnetic moment of the A+ arrives at 
fji^+ = (2.7 lb 1 lb 1.5 ± 3);Ua? . Also recent lattice calculations yield results on the magnetic dipole 
moment. The value of 1.6(3)/XAf for the A+, read off from fig. 6 of Ref. [10| at a pion mass of 
around 300 MeV lies slightly below our value of 2.32(16)/XAr. 



3.4 Electric quadrupole form factor Gei 

This is the first accurate evaluation of the quadrupole form factor, Ge2, within lattice QCD. 
The intrinsic quadrupole moment Q = rtf^GEii^), non-relativistically, is given by |^ 

Q = Jd^rY{r){3z^-r^)Y{r), (3.4) 

where xjf is the wave function of the A. A negative value, as that shown in Fig. ^ corresponds to an 
oblate deformation of the A. 

The value Ge2{0) = —0.4(14) obtained in reference is compatible with ours, but comes 
with a much larger statistical error. 



3.5 Magnetic octupole form factor Gm3 

Our signal for Gm3 is not very strong, as can be seen in Fig. This form factor has a small 
negative value, albeit with large statistical errors that, at present, do not exclude a zero value. To 
next to leading order of a chiral expansion a value of zero is expected [^. 
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Figure 2: The electric quadrupole form factor (left panel) and the magnetic octupole form factor (right 
panel). The notation is the same as in Fig.l. 
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K 


rriTt [MeV] 






Gmi(0) 


GEiiO) 


0.1554 
0.1558 
0.1562 


563(4) 
490(4) 
411(4) 


1.470(15) 
1.425(16) 
1.382(15) 


0.614(2) 
0.632(2) 
0.650(3) 


3.05(7) 
3.05(8) 
3.05(10) 


-0.6(3) 
-0.7(4) 
-0.7(4) 




135 




0.691(6) 


3.04(21) 





Table 1: Estimated values of the A mass, the charge radius, Gmi (0) and Ge2{0) for the pion masses consid- 
ered. The last row contains our results from data extrapolated to the physical value of the pion mass. 

4. Conclusions and outlook 

Our results on the electromagnetic form factors of the A baryon and their momentum depen- 
dence confirm the up to now phenomenological description, e.g. the ^^-dependence of Geo and 
Gmi are well described by a dipole ansatz. A particularly interesting result of our study is the 
non-vanishing, negative value of the electric quadrupole form factor, which is associated with an 
oblate A. 

While the statistical errors in our calculation are under control, there are several sources of 
systematic errors that we would Uke to address. The results presented here were obtained in the 
quenched approximation. A study with dynamical fermions is currently in progress. It is not clear 
how strongly the neglect of disconnected contributions affects our final results. This is a problem 
we share with most other calculations of electromagnetic form factors. Only recently, with new 
methods and faster machines, we are beginning to address the calculation of the disconnected 
contributions in a statistically controlled manner. 

While the finite volume effects are expected to be negUgible on our rather big lattice, cutoff 
effects may be significant, especially towards higher momentum transfers. 
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